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Abstract 

Searches for the lepton flavour violating decay t~ — > [i~ and the lepton 
flavour and baryon number violating decays r _ — > pfj, + [j,~ and t~ — > pfi~ have 
been carried out using proton-proton collision data, corresponding to an integrated 
luminosity of 1.0 fb , taken by the LHCb experiment at -y/s = 7TeV. No evidence 
has been found for any signal, and limits have been set at 90% confidence level 
on the branching fractions: B(r~ — > fi~ /.C pL~) < 8.0 x 10 -8 , B(t~ — > < 
3.3 x 10 -7 and B{r~ — > p[i~n~) < 4.4 x 10~ 7 . The results for the r _ — > pfi + ^ and 
r~ — > -pyT fj," decay modes represent the first direct experimental limits on these 
channels. 



Submitted to Physics Letters B 
© CERN on behalf of the LHCb collaboration, license |CT3-BY-3.0| 

^Authors are listed on the following pages. 




ii 



LHCb collaboration 



R. Aaij 40 , C. Abellan Beteta 35 -™, B. Adeva 36 , M. Adinolfi 45 , C. Adrover 6 , A. Affolder 51 , 

Z. Ajaltouni 5 , J. Albrecht 9 , F. Alessio 37 , M. Alexander 50 , S. Ali 40 , G. Alkhazov 29 , 

P. Alvarez Cartelle 36 , A.A. Alves Jr 24 ' 37 , S. Amato 2 , S. Amerio 21 , Y. Amhis 7 , L. Anderlini 17 ^ , 

J. Anderson 39 , R. Andreassen 56 , R.B. Appleby 53 , O. Aquines Gutierrez 10 , F. Archilli 18 , 

A. Artamonov 34 , M. Artuso 57 , E. Aslanides 6 , G. Auriemma 24 '™ - , S. Bachmann 11 , J.J. Back 47 , 

C. Baesso 58 , V. Balagura 30 , W. Baldini 16 , R.J. Barlow 53 , C. Barschel 37 , S. Barsuk 7 , 

W. Barter 46 , Th. Bauer 40 , A. Bay 38 , J. Beddow 50 , F. Bedeschi 22 , I. Bediaga 1 , S. Belogurov 30 , 

K. Belous 34 , I. Belyaev 30 , E. Ben-Haim 8 , G. Bencivenni 18 , S. Benson 49 , J. Benton 45 , 

A. Berezhnoy 31 , R. Bernet 39 , M.-O. Bettler 46 , M. van Beuzekom 40 , A. Bien 11 , S. Bifani 44 , 

T. Bird 53 , A. Bizzeti 17 '' 1 , P.M. Bj0rnstad 53 , T. Blake 37 , F. Blanc 38 , J. Blouw 11 , S. Blusk 57 , 

V. Bocci 24 , A. Bondar 33 , N. Bondar 29 , W. Bonivento 15 , S. Borghi 53 , A. Borgia 57 , 

T.J.V. Bowcock 51 , E. Bowen 39 , C. Bozzi 16 , T. Brambach 9 , J. van den Brand 41 , J. Bressieux 38 , 

D. Brett 53 , M. Britsch 10 , T. Britton 57 , N.H. Brook 45 , H. Brown 51 , I. Burducea 28 , A. Bursche 39 , 
G. Busetto 21 ' 9 , J. Buytaert 37 , S. Cadeddu 15 , O. Callot 7 , M. Calvi 20 '^', M. Calvo Gomez 35 '™, 

A. Camboni 35 , P. Campana 18 ' 37 , D. Campora Perez 37 , A. Carbone 14 ' c , G. Carboni 23 ' fc , 

R. Cardinale 19 '*, A. Cardini 15 , H. Carranza-Mejia 49 , L. Carson 52 , K. Carvalho Akiba 2 , 

G. Casse 51 , L. Castillo Garcia 37 , M. Cattaneo 37 , Ch. Cauet 9 , M. Charles 54 , Ph. Charpentier 37 , 

P. Chen 3 ' 38 , N. Chiapolini 39 , M. Chrzaszcz 25 , K. Ciba 37 , X. Cid Vidal 37 , G. Ciezarek 52 , 

P.E.L. Clarke 49 , M. Clemencic 37 , H.V. Cliff 46 , J. Closier 37 , C. Coca 28 , V. Coco 40 , J. Cogan 6 , 

E. Cogneras , P. Collins 37 , A. Comerma-Montells 35 , A. Contu 15 ' 37 , A. Cook 45 , M. Coombes 45 , 
S. Coquereau 8 , G. Corti 37 , B. Couturier 37 , G.A. Cowan 49 , D.C. Craik 47 , S. Cunliffe 52 , 

R. Currie 49 , C. D'Ambrosio 37 , P. David 8 , P.N.Y. David 40 , A. Davis 56 , I. De Bonis 4 , 

K. De Bruyn 40 , S. De Capua 53 , M. De Cian 39 , J.M. De Miranda 1 , L. De Paula 2 , W. De Silva 56 , 

P. De Simone 18 , D. Decamp 4 , M. Deckenhoff 9 , L. Del Buono 8 , N. Deleage 4 , D. Derkach 14 , 

0. Deschamps 5 , F. Dettori 41 , A. Di Canto 11 , F. Di Ruscio 23 ' fc , H. Dijkstra 37 , M. Dogaru 28 , 
S. Donleavy 51 , F. Dordei 11 , A. Dosil Suarez 36 , D. Dossett 47 , A. Dovbnya 42 , F. Dupertuis 38 , 
R. Dzhelyadin 34 , A. Dziurda 25 , A. Dzyuba 29 , S. Easo 48 ' 37 , U. Egede 52 , V. Egorychev 30 , 

S. Eidelman 33 , D. van Eijk 40 , S. Eisenhardt 49 , U. Eitschberger 9 , R. Ekelhof 9 , L. Eklund 50 ' 37 , 

1. El Rifai 5 , Ch. Elsasser 39 , D. Elsby 44 , A. Falabella 14 ' 6 , C. Farber 11 , G. Fardell 49 , C. Farinelli 40 , 
S. Farry 12 , V. Fave 38 , D. Ferguson 49 , V. Fernandez Albor 36 , F. Ferreira Rodrigues 1 , 

M. Ferro-Luzzi 37 , S. Filippov 32 , M. Fiore 16 , C. Fitzpatrick 37 , M. Fontana 10 , F. Fontanelli 19 '*, 
R. Forty 37 , O. Francisco 2 , M. Frank 37 , C. Frei 37 , M. Frosini 17 ' / , S. Furcas 20 , E. Furfaro 23 - fc , 
A. Gallas Torreira 36 , D. Galli 14 ' c , M. Gandelman 2 , P. Gandini 57 , Y. Gao 3 , J. Garofoli 57 , 
P. Garosi 53 , J. Garra Tico 46 , L. Garrido 35 , C. Gaspar 37 , R. Gauld 54 , E. Gersabeck 11 , 
M. Gersabeck 53 , T. Gershon 47 ' 37 , Ph. Ghez 4 , V. Gibson 46 , V.V. Gligorov 37 , C. Gobel 58 , 

D. Golubkov 30 , A. Golutvin 52 ' 30 ' 37 , A. Gomes 2 , H. Gordon 54 , M. Grabalosa Gandara 5 , 
R. Graciani Diaz 35 , L.A. Granado Cardoso 37 , E. Grauges 35 , G. Graziani 17 , A. Grecu 28 , 

E. Greening 54 , S. Gregson 46 , P. Griffith 44 , O. Griinberg 59 , B. Gui 57 , E. Gushchin 32 , 

Yu. Guz 34 ' 37 , T. Gys 37 , C. Hadjivasiliou 57 , G. Haefeli 38 , C. Haen 37 , S.C. Haines 46 , S. Hall 52 , 

T. Hampson 45 , S. Hansmann-Menzemer 11 , N. Harnew 54 , S.T. Harnew 45 , J. Harrison 53 , 

T. Hartmann 59 , J. He 37 , V. Heijne 40 , K. Hennessy 51 , P. Henrard 5 , J. A. Hernando Morata 36 , 

E. van Herwijnen 37 , E. Hicks 51 , D. Hill 54 , M. Hoballah 5 , C. Hombach 53 , P. Hopchev 4 , 

W. Hulsbergen 40 , P. Hunt 54 , T. Huse 51 , N. Hussain 54 , D. Hutchcroft 51 , D. Hynds 50 , 

V. Iakovenko 43 , M. Idzik 26 , P. Ilten 12 , R. Jacobsson 37 , A. Jaeger 11 , E. Jans 40 , P. Jaton 38 , 



iii 



F. Jing 3 , M. John 54 , D. Johnson 54 , C.R. Jones 46 , C. Joram 37 , B. Jost 37 , M. Kaballo 9 , 

S. Kandybei 42 , M. Karacson 37 , T.M. Karbach 37 , I.R. Kenyon 44 , U. Kerzel 37 , T. Ketel 41 , 

A. Keune 38 , B. Khanji 20 , O. Kochebina 7 , I. Komarov 38 , R.F. Koopman 41 , P. Koppenburg 40 , 
M. Korolev 31 , A. Kozlinskiy 40 , L. Kravchuk 32 , K. Kreplin 11 , M. Kreps 47 , G. Krocker 11 , 

P. Krokovny 33 , F. Kruse 9 , M. Kucharczyk 20 ' 25 '- 5 ; , V. Kudryavtsev 33 , T. Kvaratskheliya 30 ' 37 , 
V.N. La Thi 38 , D. Lacarrere 37 , G. Lafferty 53 , A. Lai 15 , D. Lambert 49 , R.W. Lambert 41 , 
E. Lanciotti 37 , G. Lanfranchi 18 , C. Langenbruch 37 , T. Latham 47 , C. Lazzeroni 44 , R. Le Gac 6 , 
J. van Leerdam 40 , J. -P. Lees 4 , R. Lefevre 5 , A. Leflat 31 , J. Lefrancois 7 , S. Leo 22 , O. Leroy 6 , 
T. Lesiak 25 , B. Leverington 11 , Y. Li 3 , L. Li Gioi 5 , M. Liles 51 , R. Lindner 37 , C. Linn 11 , B. Liu 3 , 

G. Liu 37 , S. Lohn 37 , I. Longstaff 50 , J.H. Lopes 2 , E. Lopez Asamar 35 , N. Lopez-March 38 , H. Lu 3 , 
D. Lucchesi 21,9 , J. Luisier 38 , H. Luo 49 , F. Machefert 7 , I.V. Machikhiliyan 4 ' 30 , F. Maciuc 28 , 

0. Maev 29 ' 37 , S. Malde 54 , G. Manca 15 ^, G. Mancinelli 6 , U. Marconi 14 , R. Marki 38 , J. Marks 11 , 
G. Martellotti 24 , A. Martens 8 , L. Martin 54 , A. Martin Sanchez 7 , M. Martinelli 40 , 

D. Martinez Santos 41 , D. Martins Tostes 2 , A. Massafferri 1 , R. Matev 37 , Z. Mathe 37 , 

C. Matteuzzi 20 , E. Maurice 6 , A. Mazurov 16 ' 32 ' 37 ' 6 , J. McCarthy 44 , A. McNab 53 , R. McNulty 12 , 

B. Meadows 56 ' 54 , F. Meier 9 , M. Meissner 11 , M. Merk 40 , D.A. Milanes 8 , M.-N. Minard 4 , 
J. Molina Rodriguez 58 , S. Monteil 5 , D. Moran 53 , P. Morawski 25 , M.J. Morello 22 ' 5 , 

R. Mountain 57 , I. Mous 40 , F. Muheim 49 , K. Midler 39 , R. Muresan 28 , B. Muryn 26 , B. Muster 38 , 

P. Naik 45 , T. Nakada 38 , R. Nandakumar 48 , I. Nasteva 1 , M. Needham 49 , N. Neufeld 37 , 

A.D. Nguyen 38 , T.D. Nguyen 38 , C. Nguyen-Mau 38 ' p , M. Nicol 7 , V. Niess 5 , R. Niet 9 , N. Nikitin 31 , 

T. Nikodem 11 , A. Nomerotski 54 , A. Novoselov 34 , A. Oblakowska-Mucha 26 , V. Obraztsov 34 , 

S. Oggero 40 , S. Ogilvy 50 , O. Okhrimenko 43 , R. 01deman 15 ' d , M. Orlandea 28 , 

J.M. Otalora Goicochea 2 , P. Owen 52 , A. Oyanguren 35 >°, B.K. Pal 57 , A. Palano 13,6 , 

M. Palutan 18 , J. Panman 37 , A. Papanestis 48 , M. Pappagallo 50 , C. Parkes 53 , C.J. Parkinson 52 , 

G. Passaleva 17 , G.D. Patel 51 , M. Patel 52 , G.N. Patrick 48 , C. Patrignani 19 '\ C. Pavel-Nicorescu 28 , 

A. Pazos Alvarez 36 , A. Pellegrino 40 , G. Penso 24 '', M. Pepe Altarelli 37 , S. Perazzini 14 ' c , 

D. L. Perego 20j , E. Perez Trigo 36 , A. Perez-Calero Yzquierdo 35 , P. Perret 5 , M. Perrin-Terrin 6 , 

G. Pessina 20 , K. Petridis 52 , A. Petrolini 19 '\ A. Phan 57 , E. Picatoste Olloqui 35 , B. Pietrzyk 4 , 
T. Pilar 47 , D. Pinci 24 , S. Playfer 49 , M. Plo Casasus 36 , F. Polci 8 , G. Polok 25 , A. Poluektov 47 ' 33 , 

E. Polycarpo 2 , A. Popov 34 , D. Popov 10 , B. Popovici 28 , C. Potterat 35 , A. Powell 54 , 

J. Prisciandaro 38 , V. Pugatch 43 , A. Puig Navarro 38 , G. Punzi 22 ' r , W. Qian 4 , J.H. Rademacker 45 , 

B. Rakotomiaramanana 38 , M.S. Rangel 2 , I. Raniuk 42 , N. Rauschmayr 37 , G. Raven 41 , 

S. Redford 54 , M.M. Reid 47 , A.C. dos Reis 1 , S. Ricciardi 48 , A. Richards 52 , K. Rinnert 51 , 
V. Rives Molina 35 , D.A. Roa Romero 5 , P. Robbe 7 , E. Rodrigues 53 , P. Rodriguez Perez 36 , 
S. Roiser 37 , V. Romanovsky 34 , A. Romero Vidal 36 , J. Rouvinet 38 , T. Ruf 37 , F. Ruffini 22 , 

H. Ruiz 35 , P. Ruiz Vails 35 ' , G. Sabatino 24 ' fc , J.J. Saborido Silva 36 , N. Sagidova 29 , P. Sail 50 , 

B. Saitta 15 ' d , V. Salustino Guimaraes 2 , C. Salzmann 39 , B. Sanmartin Sedes 36 , M. Sannino 19 '*, 
R. Santacesaria 24 , C. Santamarina Rios 36 , E. Santovetti 23 ' fc , M. Sapunov 6 , A. Sarti 18 '', 

C. Satriano 24 '™, A. Satta 23 , M. Savrie 16 - 6 , D. Savrina 30 * 31 , P. Schaack 52 , M. Schiller 41 , 

H. Schindler 37 , M. Schlupp 9 , M. Schmelling 10 , B. Schmidt 37 , O. Schneider 38 , A. Schopper 37 , 
M.-H. Schune 7 , R. Schwemmer 37 , B. Sciascia 18 , A. Sciubba 24 , M. Seco 36 , A. Semennikov 30 , 
K. Senderowska 26 , I. Sepp 52 , N. Serra 39 , J. Serrano 6 , P. Seyfert 11 , M. Shapkin 34 , 

1. Shapoval 16 - 42 , P. Shatalov 30 , Y. Shcheglov 29 , T. Shears 51 ' 37 , L. Shekhtman 33 , O. Shevchenko 42 , 
V. Shevchenko 30 , A. Shires 52 , R. Silva Coutinho 47 , T. Skwarnicki 57 , N.A. Smith 51 , E. Smith 54 ' 48 , 
M. Smith 53 , M.D. Sokoloff 56 , F.J.P. Soler 50 , F. Soomro 18 , D. Souza 45 , B. Souza De Paula 2 , 

B. Spaan 9 , A. Sparkes 49 , P. Spradlin 50 , F. Stagni 37 , S. Stahl 11 , O. Steinkamp 39 , S. Stoica 28 , 



iv 



S. Stone 57 , B. Storaci 39 , M. Straticiuc 28 , U. Straumann 39 , V.K. Subbiah 37 , L. Sun 56 , 

S. Swientek 9 , V. Syropoulos 41 , M. Szczekowski 27 , P. Szczypka 38 - 37 , T. Szumlak 26 , S. T'Jampens 4 , 

M. Teklishyn 7 , E. Teodorescu 28 , F. Teubert 37 , C. Thomas 54 , E. Thomas 37 , J. van Tilburg , 

V. Tisserand 4 , M. Tobin 38 , S. Tolk 41 , D. Tonelli 37 , S. Topp-Joergensen 54 , N. Torr 54 , 

E. Tournefier 4 ' 52 , S. Tourneur 38 , M.T. Tran 38 , M. Tresch 39 , A. Tsaregorodtsev 6 , P. Tsopelas 40 , 

N. Tuning 40 , M. Ubeda Garcia 37 , A. Ukleja 27 , D. Urner 53 , U. Uwer 11 , V. Vagnoni 14 , 

G. Valenti 14 , R. Vazquez Gomez 35 , P. Vazquez Regueiro 36 , S. Vecchi 16 , J.J. Velthuis 45 , 

M. Veltri 17 ' 9 , G. Veneziano 38 , M. Vesterinen 37 , B. Viaud 7 , D. Vieira 2 , X. Vilasis-Cardona 35 '™, 
A. Vollhardt 39 , D. Volyanskyy 10 , D. Voong 45 , A. Vorobyev 29 , V. Vorobyev 33 , C. VoB 59 , 

H. Voss 10 , R. Waldi 59 , R. Wallace 12 , S. Wandernoth 11 , J. Wang 57 , D.R. Ward 46 , N.K. Watson 44 , 
A.D. Webber 53 , D. Websdale 52 , M. Whitehead 47 , J. Wicht 37 , J. Wiechczynski 25 , D. Wiedner 11 , 
L. Wiggers 40 , G. Wilkinson 54 , M.P. Williams 47 ' 48 , M. Williams 55 , F.F. Wilson 48 , J. Wishahi 9 , 
M. Witek 25 , S.A. Wotton 46 , S. Wright 46 , S. Wu 3 , K. Wyllie 37 , Y. Xie 49 ' 37 , F. Xing 54 , Z. Xing 57 , 
Z. Yang 3 , R. Young 49 , X. Yuan 3 , O. Yushchenko 34 , M. Zangoli 14 , M. Zavertyaev 10 ' a , F. Zhang 3 , 
L. Zhang 57 , W.C. Zhang 12 , Y. Zhang 3 , A. Zhelezov 11 , A. Zhokhov 30 , L. Zhong 3 , A. Zvyagin 37 . 

1 Centro Brasileiro de Pesquisas Fisicas ( CBPF), Rio de Janeiro, Brazil 

2 Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 

3 Center for High Energy Physics, Tsinghua University, Beijing, China 
4 LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 

5 Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
6 CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
7 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

8 LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 
9 Fakultdt Physik, Technische Universitdt Dortmund, Dortmund, Germany 

10 Max-Planck-Institut fur Kernphysik (MPIK), Heidelberg, Germany 

11 Physikalisches Institut, Ruprecht-Karls-Universitdt Heidelberg, Heidelberg, Germany 

12 School of Physics, University College Dublin, Dublin, Ireland 
13 Sezione INFN di Bari, Bari, Italy 

14 Sezione INFN di Bologna, Bologna, Italy 

15 Sezione INFN di Cagliari, Cagliari, Italy 
16 Sezione INFN di Ferrara, Ferrara, Italy 
17 Sezione INFN di Firenze, Firenze, Italy 

18 Laboratori Nazionali dell 'INFN di Frascati, Frascati, Italy 

19 Sezione INFN di Genova, Genova, Italy 

20 Sezione INFN di Milano Bicocca, Milano, Italy 

21 Sezione INFN di Padova, Padova, Italy 

22 Sezione INFN di Pisa, Pisa, Italy 

23 Sezione INFN di Roma Tor Vergata, Roma, Italy 
24 Sezione INFN di Roma La Sapienza, Roma, Italy 

25 Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland 
26 AGH - University of Science and Technology, Faculty of Physics and Applied Computer Science, 
Krakow, Poland 

27 National Center for Nuclear Research (NCBJ), Warsaw, Poland 

28 Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

29 Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 

30 Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

31 Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 

32 Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 

33 Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 



V 



34 Institute for High Energy Physics (IHEP), Protvino, Russia 

35 Universitat de Barcelona, Barcelona, Spain 

36 Universidad de Santiago de Compostela, Santiago de Compostela, Spain 

37 European Organization for Nuclear Research (CERN), Geneva, Switzerland 
38 Ecole Polytechnique Federate de Lausanne (EPFL), Lausanne, Switzerland 
39 Physik-Institut, Universitat Zurich, Zurich, Switzerland 

40 Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 

41 Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The 

Netherlands 

42 NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 
^Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 
44 University of Birmingham, Birmingham, United Kingdom 
45 H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

46 Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 

47 Department of Physics, University of Warwick, Coventry, United Kingdom 
4S STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

49 School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 
50 School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 

51 Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 

52 Imperial College London, London, United Kingdom 

53 School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 
54 Department of Physics, University of Oxford, Oxford, United Kingdom 
55 Massachusetts Institute of Technology, Cambridge, MA, United States 

56 University of Cincinnati, Cincinnati, OH, United States 

57 Syracuse University, Syracuse, NY, United States 

58 Pontificia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to 2 
59 Institut fur Physik, Universitat Rostock, Rostock, Germany, associated to 11 

a P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia 

b Universita di Bari, Bari, Italy 

c Universita di Bologna, Bologna, Italy 

d Universita di Cagliari, Cagliari, Italy 

e Universita di Ferrara, Ferrara, Italy 

f Universita di Firenze, Firenze, Italy 

9 Universita di Urbino, Urbino, Italy 

h Universita di Modena e Reggio Emilia, Modena, Italy 

1 Universita di Genova, Genova, Italy 

Universita di Milano Bicocca, Milano, Italy 
Universita di Roma Tor Vergata, Roma, Italy 

1 Universita di Roma La Sapienza, Roma, Italy 
m Universita delta Basilicata, Potenza, Italy 

™ LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain 

°IFIC, Universitat de Valencia- C SIC, Valencia, Spain 

p Hanoi University of Science, Hanoi, Viet Nam 

q Universita di Padova, Padova, Italy 

r Universita di Pisa, Pisa, Italy 

s Scuola Normale Superiore, Pisa, Italy 



vi 



1 Introduction 



The observation of neutrino oscillations was the first evidence for lepton flavour violation 
(LFV). As a consequence, the introduction of mass terms for neutrinos in the Standard 
Model (SM) implies that LFV exists also in the charged sector, but with branching fractions 
smaller than ~ 10~ 40 [TJ[2] . Physics beyond the Standard Model (BSM) could significantly 
enhance these branching fractions. Many BSM theories predict enhanced LFV in r~ 
decays with respect to \x~ decays^} with branching fractions within experimental reach J3]. 
To date, no charged LFV decays such as /i~ —> e~j, \T —> e~e + e~, r~ — > £~j and 
r~ — > (with i~ = e~,fi~) have been observed j4j. Baryon number violation (BNV) 

is believed to have occurred in the early universe, although the mechanism is unknown. 
BNV in charged lepton decays automatically implies lepton number and lepton flavour 
violation, with angular momentum conservation requiring the change \A(B — L)\ = 
or 2, where B and L are the net baryon and lepton numbers. The SM and most of its 
extensions [I] require \A(B — L)\ =0. Any observation of BNV or charged LFV would 
be a clear sign for BSM physics, while a lowering of the experimental upper limits on 
branching fractions would further constrain the parameter spaces of BSM models. 

In this Letter we report on searches for the LFV decay r~ — > [i + fi~ and the LFV 
and BNV decay modes t~ — > pfi + pr and r~ — > p^i~/j,~ at LHCb pi. The inclusive r~ 
production cross-section at the LHC is relatively large, at about 80 U-b (approximately 
80% of which comes from D~ — > t~v t ), estimated using the bb and cc cross-sections 
measured by LHCb (6j[7] and the inclusive b — > r and c — > r branching fractions [8j. The 
r _ — > fi~fi + fi~ and r — > pn/i decay modes^jare of particular interest at LHCb, since muons 
provide clean signatures in the detector and the ring-imaging Cherenkov (RICH) detectors 
give excellent identification of protons. 

This Letter presents the first results on the r~ — > fi~fi + fi~ decay mode from a hadron 
collider and demonstrates an experimental sensitivity at LHCb, with data corresponding to 
an integrated luminosity of l.Ofb -1 , that approaches the current best experimental upper 
limit, from Belle, B{r~ ->• fi~V + V~) < 2- 1 x !0 -8 at 90% confidence level (CL) BaBar 
and Belle have searched for BNV r decays with \ A(B — L)\ = and \ A(B—L)\ = 2 using the 
modes r~ — > Ah~ and Ah~ (with h~ = tt~ , K~), and upper limits on branching fractions of 
order 10~ 7 were obtained |4j. BaBar has also searched for the B meson decays B° — > A+l~ } 
B~ At (both having | A (B - L)\ = 0) and B~ Al~ (|A(S - L)\ = 2), obtaining 
upper limits at 90% CL on branching fractions in the range (3.2 — 520) x 10~ 8 |10|. The 
two BNV t decays presented here, r _ — > pfi + fi~ and t~ — > pfi~fi~, have \A(B — L)\ = 
but they could have rather different BSM interpretations; they have not been studied by 
any previous experiment. 

In this analysis the LHCb data sample from 2011, corresponding to an integrated 
luminosity of l.Ofb - collected at ^fs = 7TeV, is used. Selection criteria are implemented 
for the three signal modes, r~ — > fi~ fi + fi~ , t~ — > pfi + and r~ — > p^T , and for the 
calibration and normalisation channel, which is — > <pir~ followed by — > , referred 

1 The inclusion of charge conjugate processes is implied throughout this Letter. 

2 In the following t — > refers to both the t~ — > p/j, + fi~ and r _ — > p)i~ n~ channels. 
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to in the following as D~ — > (f)(fi + . These initial, cut-based selections are designed 
to keep good efficiency for signal whilst reducing the dataset to a manageable level. To 
avoid potential bias, and p/z/z candidates with mass within ±30MeV/c 2 (~ 3<r m ) 

of the r mass are initially blinded from the analysis, where a m denotes the expected mass 
resolution. For the 3/x channel, discrimination between potential signal and background is 
performed using a three-dimensional binned distribution in two likelihood variables and the 
mass of the r candidate. One likelihood variable is based on the three-body decay topology 
and the other on muon identification. For the r — > pfifi channels, the use of the second 
likelihood function is replaced by cuts on the proton and muon particle identification (PID) 
variables. The analysis strategy and limit-setting procedure are similar to those used for 
the LHCb analyses of the B® — > fi + n~ and B° — > fi + n~ channels 11, 12 



2 Detector and triggers 



The LHCb detector [5] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r] < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon- 
strip detectors and straw drift tubes placed downstream. The combined tracking system 
has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and impact parameter resolution of 20 |im for tracks with high transverse momentum 
(Pt)- Charged hadrons are identified using two RICH detectors. Photon, electron and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers. 



The trigger 13 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage that applies a full event reconstruction. 
The hardware trigger selects muons with p-r > 1.48GeV/c. The software trigger requires a 
two-, three- or four-track secondary vertex with a high sum of the of the tracks and a 
significant displacement from the primary pp interaction vertices (PVs). At least one track 
should have > 1-7 GeV/c and impact parameter chi-squared (IP x 2 )> with respect to the 
pp collision vertex, greater than 16. The IP x 2 is defined as the difference between the \ 2 
of the PV reconstructed with and without the track under consideration. A multivariate 
algorithm is used for the identification of secondary vertices. 

with a specific 



For the simulation, pp collisions are generated using Pythia 6.4 14 



LHCb configuration [15]. Particle decays are described by EvtGen 16 in which final- 
state radiation is generated using Photos 17 . For the three signal r decay channels, the 



final-state particles are distributed according to three-body phase space. The interaction 
of the generated particles with the detector, and its response, are implemented using the 



Geant4 toolkit 1 18 1 as described in Ref. 19 
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3 Signal candidate selection 



The signal and normalisation channels have the same topology, the signature of which is a 
vertex displaced from the PV, having three tracks that are reconstructed to give a mass 
close to that of the r lepton (or D s meson for the normalisation channel). In order to 
discriminate against background, well- reconstructed and well-identified muon, pion and 
proton tracks are required, with selections on track quality criteria and a requirement 
of pt > 300 MeV/c. Furthermore, for the r — > pfifi signal and normalisation channels 
the muon and proton candidates must pass loose PID requirements and the combined 
of the three-track system is required to be greater than 4 GeV/ c. All selected tracks are 
required to have IP x 2 > 9. The fitted three-track vertex has to be of good quality, with a 
fit x 2 < 15, and the measured decay time, t, of the candidate forming the vertex has to be 
compatible with that of a heavy meson or tau lepton (ct > 100 urn). Since the Q- values 
in decays of charm mesons to r are relatively small, poorly reconstructed candidates are 
removed by a cut on the pointing angle between the momentum vector of the three-track 
system and the line joining the primary and secondary vertices. In the r~ — > fi~fi + fi~ 
channel, signal candidates with a fi + fi~ mass within ±20MeV/c 2 of the 4> meson mass are 
removed, and to eliminate irreducible background near the signal region arising from the 
decay D~ — > ^(/i + /i~7) / u~i/ M , candidates with a mass combination below 450MeV/c 2 

are also rejected (see Section [6]). Finally, to remove potential contamination from pairs of 
reconstructed tracks that arise from the same particle, same-sign muon pairs with mass 
lower than 250MeV/c 2 are removed in both the r~ — > fi + pr and r~ — > p^T channels. 
The signal regions are defined by ±20MeV/c 2 (~ 2cr m ) windows around the nominal r 
mass, but candidates within wide mass windows, of ±400MeV/c 2 for r~ — > fi~ fi + {i~ decays 
and ±250MeV/c 2 for r — > p\i\i decays, are kept to allow evaluation of the background 
contributions in the signal regions. A mass window of ±20MeV/c 2 is also used to define 
the signal region for the D~ — » (fi(fi + channel, with the fi + fi~ mass required to be 
within ±20MeV/c 2 of the <fi meson mass. 



4 Signal and background discrimination 

After the selection each r candidate is given a probability to be signal or background 
according to the values of several likelihoods. For r~ — > fi + pr three likelihoods are used: 
a three-body likelihood, A^body, a PID likelihood, A-^pid, and an invariant mass likelihood. 
The likelihood A^ody uses the properties of the reconstructed r decay to distinguish 
displaced three-body decays from A-body decays (with N > 3) and combinations of tracks 
from different vertices. Variables used include the vertex quality and its displacement from 
the PV, and the IP and fit x 2 values of the tracks. The likelihood A^pid quantifies the 
compatibility of each of the three particles with the muon hypothesis using information 
from the RICH detectors, the calorimeters and the muon stations; the value of A^pid is 
taken as the smallest one of the three muon candidates. For r — > pfifi, the use of A^pid is 
replaced by cuts on PID quantities. The invariant mass likelihood uses the reconstructed 
mass of the r candidate to help discriminate between signal and background. 
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For the A^3bod y likelihood a boosted decision tree 120] is used, with the AdaBoost 
algorithm 21 , and is implemented via the TMVA 22 toolkit. It is trained using signal and 
background samples, both from simulation, where the composition of the background is a 
mixture of bb — > fifiX and cc — > fifiX according to their relative abundance as measured 
in data. The A^pid likelihood uses a neural network, which is also trained on simulated 
events. The probability density function shapes are calibrated using the D~ — > /i~)vr~ 
control channel and J/ij) — > fi + fi~ data for the A^body and A^pid likelihoods, respectively. 
The shape of the signal mass spectrum is modelled using D~ — > 0(/i + /i~)7r~ data. The 
A^sbody response as determined using the training from the t~ — > fi~fi + [i~ samples is used 
also for the r — > pfifi analyses. 

For the A^ody an d A^pid likelihoods the binning is chosen such that the separation 
power between the background-only and signal-plus-background hypotheses is maximised, 
whilst minimising the number of bins. For the Alsbody likelihood the optimum number of 
bins is found to be six for the r~ — > ^ pT analysis and five for r — > pup, while for the 
Alpio likelihood the optimum number of bins is found to be five. The lowest bins in Al3bod y 
and AIpid do not contribute to the sensitivity and are later excluded from the analyses. 
For the r — > pfifi analysis, further cuts on the muon and proton PID hypotheses are 
used instead of AIpid and are optimised, for a 2a significance, on simulated signal events 



and data sidebands using the figure of merit from Ref. 23 , with the distributions of the 
PID variables corrected according to those observed in data. The expected shapes of the 
invariant mass spectra for the r _ — > fi~[i + fi~ and r — > p/i/i signals, with the appropriate 
selections applied, are taken from fits to the D~ — > <f)(fi + fi^)-K~ control channel in data 
as shown in Fig. [T] The signal distributions are modelled with the sum of two Gaussian 
functions with a common mean, where the narrower Gaussian contributes 70% of the total 
signal yield, while the combinatorial backgrounds are modelled with linear functions. The 
expected widths of the r signals in data are taken from simulation, scaled by the ratio 
of the widths of the D~ peaks in data and simulation. The data are divided into eight 
equally spaced bins in the ±20 MeV/ c 2 mass window around the nominal r mass. 



5 Normalisation 

To measure the signal branching fraction for the decay r~ — > fi~fi + fi~ (and similarly for 
r — > p\i\i) we normalise to the D~ — > pr)-K~ calibration channel using 



B{t — > \1 jl + jl ) 



fD a REC&SEL ,TRIG jvr 

Jt w e cal w fc cal w iv sig 



B(D- -> T-V T ) X gR^C&SEL X e TPUG >< Ncai 



= a x N sig , (1) 

where a is the overall normalisation factor and iV S ig is the number of observed signal 
events. The branching fraction B(D~ — > t~v t ) is taken from Ref. (24] . The quantity f® s 
is the fraction of r leptons that originate from D~ decays, calculated using the bb and cc 
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Figure 1: Invariant mass distribution of 0(/i + fi~ 
(b) the t —7- pfifj, selection and PID cuts. The solid (blue) lines show the overall fits, the long 
dashed (green) and short dashed (red) lines show the two Gaussian components of the signal 
and the dot dashed (black) lines show the backgrounds. 



cross-sections as measured by LHCb [6j[7] and the inclusive b — > r, c — > r, b — > D s and 
c — > D s branching fractions |8j. The corresponding expression for the r — > pfifi decay is 
identical except for the inclusion of a further term, e™/e™, to account for the effect of 
the PID cuts. 

The reconstruction and selection efficiencies, e REC&SEL ; are products of the detector 
acceptances for the particular decays, the muon identification efficiencies and the selection 
efficiencies. The combined muon identification and selection efficiency is determined from 
the yield of simulated events after the full selections have been applied. In the sample of 
simulated events, the track IPs are smeared to describe the secondary- vertex resolution of 
the data. Furthermore, the events are given weights to adjust the prompt and non-prompt 
b and c particle production fractions to the latest measurements M. The difference in 
the result if the weights are varied within their uncertainties is assigned as a systematic 
uncertainty. The ratio of efficiencies is corrected to account for the differences between data 
and simulation in efficiencies of track reconstruction, muon identification, the 0(1020) mass 
window cut in the normalisation channel and the r mass window cut, with all associated 
systematic uncertainties included. The removal of candidates in the least sensitive bins in 
the A^3bod y and A^pid classifiers is also taken into account. 

The trigger efficiency for selected candidates, e™ 6 , is evaluated from simulation while 
its systematic uncertainty is determined from the difference between trigger efficiencies of 
B~ — > J/ipK~ decays measured in data and in simulation. 

For the r — > pfifi channels the PID efficiency for selected and triggered candidates, 
e PID , is calculated using data calibration samples of J/tj) — > and A — > pn~ decays, 

with the tracks weighted to match the kinematics of the signal and calibration channels. 
A systematic uncertainty of 1% per corrected final-state track is assigned |7j, as well 
as a further 1% uncertainty to account for differences in the kinematic binning of the 
calibration samples between the analyses. 

The branching fraction of the calibration channel is determined from a combination of 
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known branching fractions using 



B{D- -+ 0(/*V)O = ^^k+K-^ ^ ^ mV) = (L33 ± °' 12) X 10 ^ ' 

(2) 

where B(<f) -> K + K ) and B(<f> ->■ are taken from (8) and £>(£>,, -)• <p(K + K )ir~ 



is taken from the BaBar amplitude analysis 25 , which considers only the (J) — > K + K~ 
resonant part of the D~ decay. This is motivated by the negligible contribution of 
non-resonant D~ — > n + events seen in our data. The yields of D~ — > 0(/i + yU~)7r~ 
candidates in data, N ca ,\, are determined from the fits to reconstructed (fi(fi + [i~)iy~ mass 
distributions, shown in Fig. [TJ The variations in the yields if the relative contributions of the 
two Gaussian components are varied in the fits are considered as systematic uncertainties. 
Table [T] gives a summary of all contributions to a; the uncertainties are taken to be 
uncorrelated. 



6 Background studies 

The background processes for the decay r~ — > \l~ p^yT consist mainly of decay chains of 
heavy mesons with three real muons in the final state or with one or two real muons in 
combination with two or one misidentified particles. These backgrounds vary smoothly 
in the mass spectra in the region of the signal channel. The most important peaking 
background channel is found to be D~ — > r](fi + / u _ 7)/i~z/ At , about 80% of which is removed 
(see Section [3]) by a cut on the dimuon mass. The small remaining background from 
this process is consistent with the smooth variation in the mass spectra of the other 



Table 1: Terms entering in the normalisation factor a for r — > fi fi + fi , r — > pfi^ fi and 
r~ — > pfi~ fi~ , and their combined statistical and systematic uncertainties. 





T — > /J, fl + fi 




t — y p/i /j, 


B(D 7 -> 0(/x+ M -)7r-) 


(1.33 ± 0.12) x 10" 5 


f?' 


0.78 ± 0.05 


B{D- r-v T ) 


0.0561 ± 0.0024 


REC&SEL / REC&SEL 
^cal / ^sig 


1.49 ±0.12 


1.35 ± 0.12 


1.36 ±0.12 


TRIG / TRIG 
^cal / ^sig 


0.753 ±0.037 


1.68 ±0.10 


2.03 ±0.13 


PID / PID 
^cal / ^sig 


n/a 


1.43 ± 0.07 


1.42 ±0.08 




48 076 ± 840 


8145 ± 180 


a 


(4.34 ±0.65) x 10~ 9 


(7.4 ± 1.2) x 10~ 8 


(9.0 ± 1.5) x 10~ 8 
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backgrounds in the mass range considered in the fit. Based on simulations, no peaking 
backgrounds are expected in the r — > pn\i analyses. 

The expected numbers of background events within the signal region, for each bin 
in -M3bod y , A^pid (for t~ — > pr p + n~) and mass, are evaluated by fitting the candidate 
mass spectra outside of the signal windows to an exponential function using an extended, 
unbinned maximum likelihood fit. The small differences obtained if the exponential curves 
are replaced by straight lines are included as systematic uncertainties. For t~ — > pT p + p~ 
the data are fitted over the mass range 1600 - 1950 MeV/c 2 , while for r ppp the fitted 
mass range is 1650 — 1900 MeV/c 2 , excluding windows around the expected signal mass of 
±30MeV/c 2 for pT p + p~ and ±20MeV/c 2 for ppp. The resulting fits to the data sidebands 
for a selection of bins for the three channels are shown in Fig. [2j 



> 



T3 

'-3 
a 



> 



6 - 



5 - 



2 4 



a 
-a 
'■5 
c 

a 
U 



(a) 



-Msbody e [0.65, 1.0] 

e [0.725,1 0] LHCb 




1600 



1700 




1800 1900 
m(ji~/^ + ju~) [MeV/c 2 ] 

10 



1700 



1800 1900 
m{pfi + jj.~) [MeV/c 2 ] 




1700 



1800 1900 
m(pjj.~jj.~) [MeV/c 2 ] 



Figure 2: Invariant mass distributions and fits to the mass sidebands in data for (a) fi + fi~ fi~ 
candidates in the four merged bins that contain the highest signal probabilities, (b) pfi + n~ 
candidates in the two merged bins with the highest signal probabilities, and (c) pfi~ fi~ candidates 
in the two merged bins with the highest signal probabilities. 
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7 Results 



Tables [2] and [3] give the expected and observed numbers of candidates for all three channels 
investigated, in each bin of the likelihood variables. No significant evidence for an excess of 
events is observed. Using the CL S method |T2] as a statistical framework, the distributions 
of observed and expected CL S values are calculated as functions of the assumed branching 
fractions, as shown in Fig. [3j The expected limits at 90% (95%) CL for the branching 
fractions are 

B{t~ -> xx~/x + /x~) < 8.3 (10.2) x 10~ 8 , 
B{j- -> pfj+fj,-) < 4.6 (5.9) x 10~ 7 , 
B{j~ -tpiTiT) < 5.4 (6.9) x 10" 7 , 

while the observed limits at 90% (95%) CL are 

B(t~ -> ,ix~/x + /x~) < 8.0 (9.8) x 10~ 8 , 
B{j- ^p/x + /x~) < 3.3 (4.3) x 10" 7 , 
B{j- ->■ pfi- lT) < 4.4 (5.7) x 10" 7 . 

All limits are given for the phase-space model of r decays. For r _ — > iT ji^ iT , the 
efficiency is found to vary by no more than 20% over the lT if mass range and by 10% 
over the li + a~ mass range. For r — > p/x/x, the efficiency varies by less than 20% over the 
dimuon mass range and less than 10% with p/x mass. 

In summary, a first limit on the lepton flavour violating decay mode r~ — > u~ li + li~ 
has been obtained at a hadron collider. The result is compatible with previous limits and 
indicates that with the additional luminosity expected from the LHC over the coming 
years, the sensitivity of LHCb will become comparable with, or exceed, those of BaBar 
and Belle. First direct upper limits have been placed on the branching fractions for two 
r decay modes that violate both baryon number and lepton flavour, r~ — > pfi + fi~ and 

T~ — > p/X~ IT . 
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Table 2: Expected background candidate yields, with their systematic uncertainties, and observed 
candidate yields within the r signal window in the different likelihood bins for the r~ — > fi~fi + fi~ 
analysis. The likelihood values for .M pid range from (most background-like) to +1 (most signal- 
like), while those for 7W3body range from —1 (most background- like) to +1 (most signal- like) . 
The lowest likelihood bins have been excluded from the analysis. 



A A 


A A 

^3body 


Expected 


Observed 




— U.4o — U.Uo 


34b.U ± D. ( 


a an 

4uy 




r\ nr nor 

0.05 - 0.35 


83.8 ± 6.6 


68 


n a o n n 

0.43 - 0.6 


nor c\ r*r 

0.35 - 0.65 


on o i o n 

30.2 ± 2.0 


o r 

35 




0.65 - 0.74 


A o i r\ o 

4.3 ± 0.8 


2 




0.74 - 1.0 


1.4 ± 0.4 


1 




_n 48 — 05 


VQ 1 _|_ Q 1 
(O.J- _l_ tJ. -L 


64 




0.05-0.35 


18.3 ± 1.5 


15 


0.6 - 0.65 


0.35-0.65 


8.6 ± 1.1 


7 




0.65-0.74 


0.4 ± 0.1 







0.74-1.0 


0.6 ± 0.2 


2 




-0.48-0.05 


45.4 ± 2.4 


51 




0.05-0.35 


11.7 ± 1.2 


6 


0.65 - 0.725 


0.35-0.65 


5.3 ± 0.8 


3 




0.65-0.74 


0.8 ± 0.2 


1 




0.74-1.0 


0.4 ± 0.1 







-0.48-0.05 


44.5 ± 2.4 


62 




0.05-0.35 


10.6 ± 1.2 


13 


0.725 - 0.86 


0.35-0.65 


7.3 ± 1.0 


7 




0.65-0.74 


1.0 ± 0.2 


2 




0.74-1.0 


0.4 ± 0.1 







-0.48-0.05 


5.9 ± 0.9 


7 




0.05-0.35 


0.7 ± 0.2 


1 


0.86 - 1.0 


0.35-0.65 


1.0 ± 0.2 


1 




0.65-0.74 


0.5 ± 0.0 







0.74-1.0 


0.4 ± 0.1 
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Table 3: Expected background candidate yields, with their systematic uncertainties, and observed 
candidate yields within the r mass window in the different likelihood bins for the r — > p/i/i 
analysis. The likelihood values for A^3body range from —1 (most background-like) to +1 (most 
signal- like) . The lowest likelihood bin has been excluded from the analysis. 





r — > pn + fi 


r — )■ pfj, /j, 


-A^3body 


Expected 


Observed 


Expected 


Observed 


-0.05-0.20 


37.9 ± 0.8 


43 


41.0 ± 0.9 


41 


0.20-0.40 


12.6 ±0.5 


8 


11.0 ±0.5 


13 


0.40 - 0.70 


6.76 ± 0.37 


6 


7.64 ± 0.39 


10 


0.70-1.00 


0.96 ± 0.14 





0.49 ±0.12 







20 40 60 20 40 60 

BR(t~— > p n + ff) [x 10" 8 ] BR(t~— > p ff IT) [x 10 -8 ] 



Figure 3: Distribution of CL S values as functions of the assumed branching fractions, under the 
hypothesis to observe background events only, for (a) t~ — > fi~ fi + fi~ , (b) t~ — > pp, + fjT and (c) 
r _ — > pfj,~pi~. The dashed lines indicate the expected curves and the solid lines the observed 
ones. The light (yellow) and dark (green) bands cover the regions of 68% and 95% confidence for 
the expected limits. 
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